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ABSTRACT 

Electron scale solar wind turbulence has attracted great interest in recent years. Clear evidences 
have been given from the Cluster data that turbulence is not fully dissipated near the proton scale but 
continues cascading down to the electron scales. However, the scaling of the energy spectra as well 
as the nature of the plasma modes involved at those small scales are still not fully determined. Here 
we survey 10 years of the Cluster search-coil magnetometer (SCM) waveforms measured in the solar 
wind and perform a statistical study of the magnetic energy spectra in the frequency range [1, 180]Hz. 
We show that a large fraction of the spectra exhibit clear breakpoints near the electron gyroscale pe, 
followed by steeper power-law like spectra. We show that the scaling below the electron breakpoint 
cannot be determined unambiguously due to instrumental limitations that will be discussed in detail. 
We compare our results to those reported in other studies and discuss their implication on the physical 
mechanisms and the theoretical modeling of energy dissipation in the SW. 
Subject headings: turbulence, solar wind, heating, whistler, KAW 



1. INTRODUCTION 

The Solar Wind (SW) is certainly the astrophysical 
plasma in which significant progress in understanding 
plasma turbulence has been most achieved. The main 
reason for that is the availability of high quality data 
from several space missions, from the earlier spacecraft 
Voyager 1 and 2 to more recent ones such as Cluster, 
Wind or Stereo. While earlier spacecraft data have al- 
lowed for addressing turbulence at scales larger than the 
ion scale (typically pi ~ 100 km, which corresponds to an 
observed frequency in the spacecraft frame of fsc ^ 0.5 
Hz using the Taylor assumption), i.e. the so-called in- 
ertial range, the high time resolution Cluster data have 
recently opened a new chapter in turbulence studies that 
focuses on electron scales {L ~ p^. ^ 1 km) (Sahraoui. 
eFal. 2009; Kiyani etfal. ' 2009 ; Ale xandrova et al. |piMf 
^hraoui ct al.,,2010a| ). These observations have driven, 
intensive research work on electron scale turbulence 



both theoretically a nd numeric ally ( Chani 



Meyrand fc Galtier 2010; Pode sta et al. 
ale & Burgess"2011 



2012| [Meyrand &: Galtier 2 0l2| : 



Howes et aL] 2011a 

Determining the nature 



et aL||2011 
0[|Gampore- 
[Sahr aoui et al 



and properties (e.g., scaling, anisotropy) of the turbu- 
lence at small scales is indeed crucial to understanding 
the problems of energy dissipation and heating, particle 
acceleration, and magnetic reco nnection in spa ce and as- 
trophysical plasmas (Schekochihin et al. 2009). Thanks 
to the high time resolution Search-Coil Magnetometer 
(SCM), which provides waveforms of the magnetic field 
fluctuations up to 180Hz, it has become clear that SW 
turbulence does cascade below the ion scale pi down to 
the electron scale pe where dissipation becomes impor- 
tant and the spectra steepen significantly. However, the 
underlying physics is still very controversial, owing to 
the fact that the available observations are very few and 
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the theoretical and numerical work has b een extended 
to those small scales o nly very recently ([Chang et al. 



20111 Camporea le & B urgess 



Sahraoui et al.[[2012( |Meyranc 



2011 



Howes et al. 2011a| 



& Galtier||2012|). 



ii'rom the observational point ot view only a handful of 
observations exist on electron scales SW turbulence and 
all do not agree. The first point of controversy is the scal- 
ing o£_Hie_jna£nBtic_ener spectra down to and below 
Pe- Sahraoui et al. (2009) have first reported a power-law 
cascade /"^ '^ down to fp^ , where a clear spectral break is 
observed and followed by another power-law like shown 
to be close to /sc~^ (/p^ is the frequency in the spacecraft 
frame corresponding to the electron gyroscale when the 
Taylor frozen- in-fiow assumption is usecQ. It was em- 
phasized in that reference that the fit below fp^ is subject 
to caution owing to the limited extension of the elec- 
tron dissipation range, imposed by the sensitivity of the 
SCM (see section [2| belo w). Similar power-law-like spec- 



et al. 



tra ha v e been reported in Sahraoui et al.l(l2010al);l Kiyani 
( 2009 1 . On the other hand pVlcxandrova et al 



1 2ro| ,llsmg STAFF-Spectrum Analyzer (STAFF^SX^ 
data, have reported exponential scaling ^ expi —\/kpe) 
in the scale r ange fcpe ~ [0.1,1]. More recently Alexan- 



drova et al. (2012) proposed a different model where a 

scaling kj_~ ^ exp{—k±pe) was suggested to fit better the 
data in the scale range kp^ ^ [0.03,3]. We will refer to 
the above three models respectively by the double-power- 
law, the exponential and the hybrid. Besides these three 
models we introduce below a new one, the asymptotic- 
double-power-law. We show that the exponential model 
does not fit well the data as does the double-power-law 
model. We show however that, due to Cluster instru- 



^ It is commonly referred to as the Doppler-shifted electron gyro- 
radius, although it is not a very appropriate definition considering 
that Doppler shift applies to frequencies and not to spatial scales. 
Here we introduce the terminology of the "Taylor-shifted scale" . 



2 



Sahraoui et al. 



mental limitations that we shall discuss in detail, this 
latter may not be always possible to distinguish within 
the data from the two remaining models, the main dif- 
ference between them being reduced to the presence or 
not of the spectral break near p^. 

From the theoretical point of view, very recent numer- 
ical simulations have tried to tackle the problem of cas- 
cade and dissip ation at electrons scales. |Camporeale fc 



Burgess] ( |2011[ ) performed 2D full Particle-In-Cell (PIC) 
simulations of decaying electromagnetic fluctuations near 
and below the electron inertial length d^. They showed 
that the magnetic energy spectra steepen from kj^'^ to 
kj^'^ with a clear spectral break at fcpe ~ 1. 3D PIC 
simulations of whistler-driven turbulence showed very 
similar res ults with less steep spectra, k±~^'^, at scales 
k_Lde > 1 dChang et aL]|2011| ). GyroKinetic (GK) sim- 
ulations of strong Kmetic Alfven Wave (KAW) turbu- 
lence, which self-consistently contain kinetic damping of 
the low frequency plasma modes (w << uid) but the cy- 
clotron resonance, s howed clear power- law cascade fcj^'^ 
down to fc_LPe ~ 1 (Howes et al.|20lTa). The GK theory 



predicts magnetic energy spectra with a sc aling k , ^^^^ 
for th e (entropy) cascade b elow pe (^ chckoc hihin et al. 
2009|. On the other hand, IMeyrand & Galtierl (|201U|) 



showed that incompressible Electron MMD (EMMD) tur 
bulence predicts that the magnetic energy spectra should 
follow a power law k~^^/^ at scales smaller than de- That 
steepening was proposed to explain the observe d spectra 



below the electron scale reported in the SW (Sahraoui 
et al.||2009| |2010a[ ). Although this fluid model is non- 
dissipative and does not consider any damping of the 
turbulence via kinetic effects, which are important in the 
dispersive and the dissipation ranges, it has neverthe- 
less the merit of proving that fully nonlinear dynamics 
and dispersive effects in EMHD leads to a steepening 
of the energy spectra near de- Thus, all existing the- 
oretical and numerical predictions support the double- 
power-law model to reproduce the scaling near and be- 
low the electron g yrosc ale or inertial length. In a recent 
paper by plowes et al. ( |2011c[ ) it was pointed out that 
for Pi « 1 the kinetic damping of the KAW turbulence 
becomes strong, which results in a cut-off of the spectra 
at moderate wavenumber, typicall y k±Pi ^ 10. This is 
inconsistent with SW observati ons ( Sahraoui et al.|2(j09 



2010a[[Ale^androva et al.|2009p . To overcome the failure 



of the model to reproduce SW observations it was sug- 
gested in that paper that nonlocal effects, due to large 
scale shear flows, play a role in sustaining the energy 
cascade down to smaller scales. When such nonlocal ef- 
fects are considered it was indeed shown that the energy 
cascade in low (3i can reach scales k±pe ~ 1, where the 
spectra curve and show an exponential-like behaviour. 
For Pi }^ 1 relevant to the SW both the local and the 
nonlocal models give the same power-law scal ing of the 
m agnetic e nergy spectra down to k±pe ~ 1 (|Howes et 



ar]|2011a|c ) , which may indicate that nonlocal effects do 



not play an important role in this case . We will return 
to this discussion below. 

Here we report a large statistical survey of SW turbu- 
lence using the Cluster SCM data sampled at 450Hz. We 
focus particularly on the scaling of the magnetic energy 
spectra up to and above /p^, the frequency correspond- 



ing to Taylor-shifting the electron gyroscale pe- In sec- 
tion [2] we discuss several caveats that need to be handled 
carefully in order to avoid any misinterpretation of the 
observations. In section [3] we explain the approach we 
used to select the Cluster data in the SW. We particu- 
larly target intervals when the Cluster fleet is located in 
the SW, excluding thus data from the electron and ion 
foreshock regions, and when the Signal-to-Noise-Ratio 
(SNR) of the magnetic fluctuations is high. In section |5] 
we show the main results of the study and discuss them 
in light of existing theoretical and numerical predictions 
and earlier work related to the subject. 

2. EXPERIMENTAL CAVEATS 

2.1. Effect of the limited sensitivity of the SCM 

The Cluster SCM data are so far the only available 
data allowing one to probe into electron scales of SW 
turbulence with relatively high SNR. The SNR (in dB) is 
defined as a fu nction of the measured frequency onboard 
the spacecraft ( Sahraoui et al.||2"010c ) 



SNRifs,) = lOLog, 



&B\f,e) 



SB, 



\fsc) 



where 5B and SBsens. are respectively the amplitude of 
the magnetic fluctuation and the level of the sensitiv- 
ity floor of the Cluster SCM at the frequency fsc- Cur- 
rent space missions (e.g. Stereo, Wind, Themis) and the 
planned ones (e.g. MMS) either do not have SCM or do 
have but with limited sensitivity. The sensitivity of the 
SCM is essentially (but not exclusively) limited by the 
length of the magnetic sensor, which is generally con- 
strained by the payload on the spacecraft. For instance 
the SCMs onboard Cluster, Themis and MMS have sen- 
sors with respective lengths of 27cm, 15cm and 10cm 
which yields respective sensitivity levels of 0.43, 0.70, and 
2.35pT/\/ i/z@10Hz. Despite this relatively high sensi- 
tivity of the Cluster SCM compared to all other space 
missions, it does not allow us all the time to measure 
the SW magnetic fluctuations above ~ 30Hz due to their 



ver y low a mplitude (see Figs. 7 and 8 in , Sahraoui et 
al. (2010c)). Indeed, when the measured SNR is excep- 
tionally high (typically SNR> 30) one should suspect the 
crossin gs of other regions su ch as the ion or electron fore- 
shock ( Sahraoui et al. 2009 ) or boundaries such as inter- 
planetary shocks or CMFs (e.g . the event of 2004-01-22 
studied in |Alexandrova et al."] ( |2009 , 2012)). We recall 
that Cluster, i'hemis and MMS are magnetospheric mis- 
sions and have been designed to address essentially the 
physics of different regions of the Earth's magnetosphere 
and the magnetosheath where the magnetic fluctuations 
have generally higher amplitudes than in the SW. 

It follows from this discussion that in order to ad- 
dress properly the problem of energy cascade and dis- 
sipation in the SW at electron scales, typically at fre- 
quencies fsc ^ 30Hz, it is important to select SW data 
when the fluctuations amplitude is very high compared to 
the sensitivity floor of th e Cluster SCM. Because o f this 
limited SNR in the SW [Alexandrova et al.| ( p009| pro- 
posed to subtract the sensitivity level ot the instrument 
from the measured power spectra. Without discussing 
in detail the foundation of a such procedure from the 
signal processing point of view (in particular since the 
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noise of the instrument is not known accurately, and is 
generally estimated from measurements in the magneto- 
spheric lobes), we only underline how affected would be 
the scaling of the resulting "denoised" spectra. Indeed, 
when the SNR becomes very low at a given frequency 
/o, i.e. SB^f > /o) = (1 + e{f))dBsens\f > fo) with 
e(/) << 1, subtracting the sensitivity of the instrument 
from the actual spectrum yields a spectrum which falls 



to zero: SB^ — SB. 



eSB, 



It is not surpris- 



ing that on the Log-Log scale such a vanishing spectrum 
will be seen as exponentially falling-off for all frequencies 
/>/o. 

2.2. Effect of the limited frequency resolution of 
STAFF-Spectrum Analyzer 

The STAFF instrument onboard Cluster has a single 
tri-axial Search-Coil Magnetometer (SCM), but has two 
sub-experiments: STAFF-SC measures magnetic field 
fluctuations (i.e. waveforms) with two possible sampling 
time: 25 samples/sec (Normal Mode, NM) and 450 sam- 
ples/sec (Burst Mode, BM) with however a low-pass fil- 
ter cut-off at ISOHz, and STAFF-SA computes onboard 
the spacecraft the 5x5 spectral matrix of the three 
components of the magnetic field and the two compo- 
nents of the ele ctric field fluctuations measured by EFW 
( | Cornilleau eTal . . 2003) . STAFF-SA transmits the com- 
puted correlation matrix on the central frequencies fo 
of 27 channels logarithmically spaced in the frequency 
range [SHz, 4kIIz] whose bandwidth is proportional to 
fo, Af = 13%/o. In BM the interval reduces to ap- 
proximately [70Hz,4kIIz]. It is important to emphasize 
here that due to the limited SNR in the SW, which 
generally does not allow us to address frequencies (in 
the spacecraft frame) higher than ~ lOOHz and consid- 
ering that the Taylor-shifted electron gyroscale is typi- 
cally fp^ - 70H^ , STAFF-SA provides thus informa- 
tion only at a few frequencies in the interval [8, lOOjHz 
which are ~[9, 11, 14, 18, 22, 28, 35, 44, 55, 70, 88]Hz. This 
limited frequency resolution of the STAFF-SA spectra in 
the scale range between proton and electron gyroscales 
may result either in missing the spectral breaks or in not 
localizing them properly. In both cases determination of 
the scaling of the spectra at electron scales may be not 
very accurate. To show such a difficulty to locate the 
spectral break we plotted in Figjl] the same spectrum as 
in Figjsjjc) (measured by STAFF-SC) as it would have 
been measured by STAFF-SA. We can see that while 
the spectral break can be observed easily on Figj8](c) (at 
/ ^ SOHz) it becomes more difficult to identify in Fig{l] 
The limited number of available points (^ 10) would also 
make any fit subjet to high uncertainties (i.e. large error 
bars) . 

This effect may be particularly relevant to explain the 
very low correlation (C = 0.03) between the spec tral 
breaks and pe reported in (Alexandrova et al. 2012) as 
compared to the one shown in b'ig 10 (C = 0.72), which is 



obtained using the STAFF-SC BJVTwaveforms. STAFF- 
SA is certainly useful to provide information on the over- 
all shape of the spectra at high frequency. However, it 

^ For SW speed of ~ 600 km/s. Faster SW yields higher values 
of /pj, that can be larger than lOOHz. In this case using STAFF- 
SA spectra may prove to be useful to study the electron dissipation 
range forming above fp^ . 
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Fig. 1. — The spectrum of FigjSfc) as it would have been mea- 
sured by STAFF-SA with the available sampling frequencies. 

may not suffice to address specific details such as the ac- 
tual scaling of the spectra at electron scales. When avail- 
able, BM STAFF-SC data is certainly more appropriate 
to address this problem. The available range of frequen- 
cies ([1, 180] Hz) is generally sufficient to address electron 
scale physics. In cases when electron scales would corre- 
spond to higher frequencies in the spacecraft frame (e.g. 
fast SW, or quasi-parallel whistler-like fluctuations) not 
accessible to STAFF-SC one then may use STAFF-SA 
spectra to probe into electron scales (provided that the 
SNR remains sufficiently high at those frequencies). 

2.3. Effect of averaging the power spectra 

Another effect that may also influence the determina- 
tion of the scaling of the turbulent spectra at electron 
scales is the time (or frequency) averaging of the spec- 
tra, which is usually done to smooth the spectra under 
the assumption of time stationarity of the fluctuations. 
The approach consists generally in dividing a given time 
series into several shorter intervals of time for which a 
power spectrum is computed. Then the resulting spectra 
are averaged to reduce the fluctuations noise in partic- 
ular at high frequencies. However, considering that the 
turbulent fluctuations are not stationary strictly speak- 
ing, averaging the spectra may result in missing some of 
the information contained in the individual spectra. This 
is particularly true regarding the presence of a spectral 
break at electron scales. For instance, let us assume that 
each given subinterval of time yields a power spectrum 
with a clear breakpoint at a given frequency in the 
spacecraft reference frame. If the value of the frequency 
fh changes from one subinterval of time to the other, 
due for instance to fluctuations either in the mean SW 
parameters (e.g. l^wj T'ii ^ej K4) that enter in estimat- 
ing the scale pe or the frequency fp^ corresponding to 
it and if the individual spectra have comparable power, 
then averaging them will result in smoothing the spectral 
breaks yielding a curved-like spectrum in the frequency 
range corresponding to all the breaks. To show this we 
plotted in Fig. [2] (top panel) three spectra computed 
over lOsec within a window of lOmn. The three spec- 
tra show clear breakpoints approximately at 20Hz, 30Hz 
and 45Hz. When the three spectra were averaged (bot- 
tom panel) they yield a spectrum without any clear break 
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Fig. 2. — Three magnetic energy spectra measured in the SW 
showing breakpoints at different frequencies indicated by the ver- 
tical lines (top). The mean spectrum shows no clear break and 
curves at those frequencies (bottom). The horizontal curve (bot- 
tom panel) shows the compensated spectrum (f)/ P(f) where 
P(/) is the used fitting function. 

in the frequency range [20, 45] Hz. The spectrum at those 
frequencies curves and can be well fit by exponential-like 
models. 

2.4. Effect of instrumental interferences at high 
frequency 

The last instrumental problem that can affect the 
proper determination of the scaling the magnetic energy 
spectra in the SW is the presence of spectral interferences 
observed at high frequency of STAFF-SC spectra. These 
interferences may be caused by the Digital Wave Pro- 
cessor (DWP) onboard Cluster. The DWP instrument 
is designed to synchronize the wave instruments of Clus- 
ter by emitting a signal allowing communication between 
them. This signal might be at the origin of the interfer- 
ences with the STAFF magnetometers (similar interfer- 
ences are also observed on FGM spectra). Examples of 
such interferences are shown in Fig. [3j 

Although the spikes caused by such interferences are 
few, their relatively high amplitudes may nevertheless 
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Fig. 3. — An example of high frequency spikes observed on the 
spectra of Bz measured by the STAFF instrument which may be 
caused by interferences coming from the DWP instrument onboard 
spacecraft 1 (top) and spacecraft 2 (bottom). The corresponding 
smoothed spectra by a mean (white) and the median (red) func- 
tions are shown. The dashed green line is the in-flight mean sensi- 
tivity floor of the STAFF instrument. 

affect the fit of the turbulent spectra at electron scales. 
Indeed, when smoothed the intense spikes lead to over 
estimating the power at the highest part of the spectra, 
which then produces shallower spectra as shown in Fig. 
[I (top). The fi gure compares smoothings by a mean and 
a median functions on windows of size logarithmically 
increasing with frequency. The effect of interferences be- 
comes more important for spectra with low SNR, which 
emphasizes the need to use the highest available SNR in 
the SW. 

From Fig[3]we can also see that the four Cluster satel- 
lites are not affected equally by the interferences. A more 
complete study (not shown here) revealed that spacecraft 
1, and particularly its component, is most affected. 
Therefore, considering that most of the studies of the 
electron scales in the SW are based on single spacecraft 
data we recommend to use data from the spacecraft 2, 3 
or 4. Depending on the level of the SNR, the number and 
amplitudes of the spikes that would appear on the spec- 
tra, the corresponding data may need to be pre-processed 
by removing the spikes prior to any analysis. 

Finally, we note that these interferences would also 
affect the spectra computed by the STAFF-SA instru- 
ment onboard the spacecraft. Owing to the fact that the 
STAFF-SA spectra are averaged over a frequency band 
centered at each of the 27 frequencies, the contribution 
of the spikes may be accounted for in the transmitted 
energy density. Since the original waveforms are lost, it 
is not possible to quantify the effect of the interferences 
on the STAFF-SA spectra. 
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Fig. 4. — A typical example of the analyzed events, (a) is the 
spectrogram of field-aligned electron energy and (b) is the electric 
field from PEACE (SCf , no data available on SC2) and WHISPER 
(SC2) showing that the spacecraft are in the free SW from 04h05 
to 05h00 (no forcshock electrons). 

3. DATA SELECTION 

In the present study we surveyed all STAFF-SC BM 
data measured in the SW by the Cluster 3 space- 
craft from 2000 to 2011. We used AMDA dat a base 
(http://cdpp-amda.cesr.fr/DDHTML/index.htmip to se- 
lect intervals ot time when the Cluster spacecraft was 
located in the SW. AMDA provides useful routines to 
search for data according to desired physical parameters 
(e.g., high or low SW speeds, plasma (3). The obtained 
list of events is then checked event by event with respect 
to WHISPER and PEACE data to remove any interval 
containing foreshock electrons. An example of the tar- 
geted SW intervals is shown in Fig|4j 

To eliminate data from the ion foreshock region (IFS), 
or from other boundaries such as Interplanetary Shocks 
(IPS) or CMEs requires additionally to use data from 
the Cluster Ion Spectrometer (CIS). The free SW is in- 
deed generally characterized by a beam of ions with an 
energy around 1 — 2 KeV. Any enlargement of the beam 
toward higher energies generally indicates a heating of 
ions across boundaries such as shocks. The CIS spec- 
trograms available at the ESA Cluster Active Archive 
(CAA) web site (http://caa.estec.esa.int/caa/home.xmlj) 
allow one to identify the ll'Ss and IPSs. However, the 
spectrograms cannot guarantee the absence of any re- 
flected ions from the shock, which can only evidenced 
by a careful analysis of the Ion Distributions Functions 
(IDFs). In the following, to eliminate the significant 
EFSs, IFSs and IPSs we rely mostly on the examina- 
tion of the Whisper, PEACE and CIS spectrograms. If 
afterwards a doubt still remains we then analyze the cor- 
responding IDFs provided at the CAA web site to check 
the presence of reflected energetic ions from the fore- 
shoc]<[^ Despite these numerous checks of the data one 
still needs to assume that if any residual reflected ions 
from the foreshock would still exist in our data (not ob- 
servable on the spectrograms or in the analyzed IDFs) 
they do not have significant impact on the electron scale 
physics, which is the main purpose of the present paper. 
Indeed, such ions when they exist would affect ion scales 
~ Pi (rather than electron scales ~ pe) through, for in- 
stance, ion plasma instabilities that might be generated 

^ The IDFs given in the ESA/CAA website may need a care- 
ful post-treatment to obtain clean IDFs in the SW. However, do- 
ing such detailed studies of the IDFs is hardly feasible for all our 
events (Mazelle ,2012J . 
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Fig. 5. — The distribution of the Signal- to-Noise- Ratio (defined 
in the text) of the Bz magnetic energy spectra at / = 30Hz. 

in the SW plasma. 

After obtaining the list of events that survive the dif- 
ferent criteria given above, we focus only on intervals 
where magnetic fluctuations have a high SNR. To do 
that we computed the magnetic energy spectra over ap- 
proximately 10s for all intervals and estimated the value 
of the SNR at the frequency f^c ^ 30Hz. With a 
SW speed Vsw ~ 500km/s and an electron gyroradius 
Pe ~ 1km, the frequency fsc = 30Hz corresponds roughly 
to the scale kpe ~ 0.4 using the Taylor assumption (other 
histograms of SNR were computed at neighboring fre- 
quencies and found to be very similar to the one shown 
here). The obtained histogram (Fig. [5| shows that at 
fsc — 30Hz most of the turbulent spectra reach the 
sensitivity floor of the Cluster SCM, which means that 
those events cannot be used to address the electron scale 
physics. To achieve this goal we need to select inter- 
vals with the highest possible SNR. However, as can be 
seen on Fig. [5] the higher is the SNR the rare are the 
events. Here we present only the results for events having 
SNR> 10, and a number of 620 was found. These events 
represented about 15% of the total events (an "event" 
refers to a total power spectrum computed over nearly 
10s period of time). 

4. FITTING MODELS 

As we discussed above, the theoretical and numerical 
studies of SW turbulence dissipation at electron scale are 
very recent. The absence of firm theoretical predictions 
on this problem and the limitations of the available data 
make it difficult to draw firm conclusions on the actual 
processes of dissipation, or even on the actual scaling of 
the energy spectra at electron scales. To illustrate that 
we plotted in Fig|6]an example of a magnetic energy spec- 
trum measured in the SW by the STAFF-SC experiment 
(black curve). We also over plotted different models to 
fit the data: 

• Double-power-law model 

P{f) = AJ-"^ [1 - H{,f - A)] + A,f-"-H{f - /,) 

• Exponential model 

Pif) = Aexpi-a^fJTb) 

• Hybrid model 

P{f) = Ar^^exp{-f/h) 
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Fig. 6. — Example of a magnetic energy spectrum of B mea- 
sured in the SW by the Cluster 2/STAFF-SC instrument (black). 
The dashed line indicates the estimated (in-flight) sensitivity floor 
of the instrument. Examples of fitting functions are plotted for 
comparison: double-power-law (red), asymptotic-double-power- 
law (green), hybrid (blue) and exponential (cyan). The horizontal 
curves are compensated spectra B"^ (f) / P{f), where P{f) is the 
corresponding fitting function. 

• Asymptotic-double-power-law model 



P{f) = A 



i + {flhf 



Above H{f) is the Heaviside 
ai — 2.6 and a2 = 4.2, and a 



function, ^ 25Hz, 
= 15. The exponential 



model has been proposed in Alexandrova et al.| (120091 
to fit the magnetic energy spectra in the scale range 
[0.1,l]fcpfQ The hybrid model is reported in ^Alexan-, 



drova et al. 1|2012 ) and combines the two previous ones 



where the power-law dominates the lowest part of the 
spectrum while the exponential function starts acting as 
frequencies approach yielding a curved-like spectrum. 
The last model is asymptotically similar to the double- 
power-law model: for / << P(/) ~ Af~"^ and for 
f»fb P{f) ~ /-"^-' Z-^-'. However, when f ^ ft 
this model is very close to the hybrid model and results 
in a curved-like spectrum. 

One should stress here that these functions are chosen 
for pure mathematical reasons and they do not neces- 
sarily reflect actual physical mechanisms of dissipation, 
although there may exist a possible connection between 
a given scaling with a particular dissipation process. For 
instance hydrodynamic turbulence is known to yield an 
exponential dissipation, while asymptotic-double-power- 
law models have been reported in fusion plasma turbu- 
lence wliBre_th£j]mite^ a key 
role ( Hasegawa et al.|[i978 Gurcan et al. 



20091. 



The comparison between the hts in I'lgjb] shows indeed 
that the double-power-law model fits better the data as 
it captures the observed spectral break. The exponential 
model is shown not to reproduce correctl y the scaling 
in the range kpe ~ [0-1,1] as reported in Alexandrova 
(2009). This model will not be considered in the 



et al. 
rest of the paper 



However, as expected, the hybrid and 



To change from frequency to wavenumber one can replace ///f, 
by kpe using the Taylor hypothesis and assuming fp^ ^ /(,. 
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Fig. 7. — The histograms of different mean SW parameters: /Si, 
ion cyclotron frequency fd, the solar wind speed Vfiow and the 
ratio of ion to electron temperatures. 

the asymptotic-double-power-law fail to fit the spectral 
break and the neighboring frequencies but they repro- 
duce the overall shape of the spectrum, in particular the 
asymptotic-double-power-law. Nevertheless, these two 
models are very close to each other, significant differences 
appear indeed only at high frequency where the noise 
floor of the instrument is reached. From this one con- 
cludes that, given that the highest SNR values available 
from the Cluster/SCM data in the SW are at best similar 
to that of Fig. [6] it is unreahstic to hope distinguishing 
between these two models in SW observations. How- 
ever, the double-power-law model can be distinguished 
from the three other (curved) models provided that the 
spectral break can be evidenced properly. As discussed 
above to detect the spectral break at electron scales re- 
quires particular caution as its presence may easily be 
smoothed by several effects. 

5. RESULTS 

For the obtained list of events we computed the histro- 
grams of the corresponding mean plasmas parameters, 
which are given in Fig. [7] We found in particular that 
most of the events have Ti > and /3i ^ 1. For these 
events we computed the corresponding power spectra of 
the magnetic fluctuations over 10s. This time length has 
been chosen to allow for covering the frequency range 
~[1, 180]Hz. To study lower frequencies longer time in- 
tervals are required, but these frequencies will not be 
addressed here. 

Some examples of the obtained spectra are plotted on 
Fig. [8j They show that the magnetic energy cascades be- 
low the ion scale and reaches the electron gyroscale where 
a clear break is evide nced, in agreement with earh er ob- 
servations reported in Sahraoui et_al. ( 2009 |2010a I . Sim- 
ilarly to Fig. 6j we plotted on I'igjSj^a) the tnree previous 
fitting models to check which one fits better the spec- 
tra. As can be seen on the compensated spectra plotted 
on the top of the figure, the double-power-law is found 
to fit better the observed spectra at the scales near and 
below the spectral break. The presence of a clear break 
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Fig. 8. — Examples of the magnetic energy spectra measured in 
the free SW showing clear spectral break at electron scales (vertical 
dashed lines). A similar description as in Fig[6] applies. 

near the electron scales invalidates the argument that the 
break at electron scales is essentially caused by the fore- 
shock electr ons present in the data analyzed in |Sahraoui 



et al. (2009). Indeed, the presence of foreshock electrons 
may generate a "bump" in t he magnetic spe ctra at that 
scale as shown in Fig. 3(a) in Sahraoui et al. ( 2009) , but 
it is not clear as to how they can generate a spectral 
break followed by a steep pow e r-law spectrum as shown 



in Fig. 3(b) of [Sahraoui et al. (2009). In any case, the 
results shown here, which were obtained from SW data 
that contained no foreshock electrons clearly prove that 
the spectral break cannot be attributed primarily to fore- 
shock electron. It most likely results from the nonlinear 
dynamics of the plasma itself. This inde ed has been ob- 
served in 2D and 3D PIC simula tions (Camporeale & 
|Burgess||2011| | Chang et al. 2(311), predicted by exist- 
mg theories on SW turbulence! S cnekochihin et al.|2009 



[Meyrand fc Galtier 2010| and o bserved also m niagne 



tosheath turbulence ( Huang et_a l. 2013). The di s agree^ 
ment with the observations of Xlexandrova et al. ( 2009 ) 
can be explained by the lack oi universality of turbulence 
at electron scales, or by the caveats discu ssed in section[2} 
Note however that the recent study by Al exandrova et aL 
j2012) showed also that 30% (out of 100) of the analyzed 
Sl'Al'F-SA spectra agreed with the double-power-law (or 
the break) model discussed here. 

Following the arguments given here we performed sim- 
ilar double-power-law fits below and above the spectral 
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Fig. 9. — The histograms of the slopes resulting from a double- 
power-law fits of the scales below (left) and above (right) the fre- 
quency /i, corresponding to the spectral break near the electron 
gyroscale. 



break to all events found to have breaks (some events, 
mostly those corresponding to fast SW, do not show 
breaks in the frequency range accessible to STAFF-SC 
The histograms of the obtained slopes are given in Fig. 
We can see that the slopes of the spectra in the dispersive 
range (i.e. [/p., /pj) cover the range ~ [—2.5, —3.1] with 
a peak at ~ —2.8. These values are close to those pre- 
dicted b y existing theoretical and nu merical simulations 
(jCarnpo rcalc & Burgess 2011; Chan g et al.||2011 Howes 



et ai.|2011a( |Boldyrev & Per ez |2012i) and to the observa- 
tions reported in | Alexandrova et al. (12012). The spectra 



below fp^ are however steeper with slopes distributed in 
the range ~ [—3.5, —5.5] and a peak at ~ —4, which 
are in general agreement with t he predictions from the 
2D and the 3D PIC sim ulations (Camporeale & Burgess 
2011 [Chang et al. 201 1|) and with the r esults ot the break 



model reported in Alexandrova et al. (2012). This wide 
distribution suggests the lack ot universality of turbu- 
lence at electron scales as compared to that in the iner- 
tial range. A similar study has been performed on mag- 
netosheath turbulence where steeper spectra above fp^ 
were observed whose slopes distrib uted in the interva l 

[—4, —7.5] and a peak near — 5.5 ( [Huang et aL]|2013[ ). 
We will discuss this point further below. 

Another important question that we can try to answer 
here is the actual physical scale, pe or de, that would 
correspond to the spectral breaks observed within the 
data. This question relates to the appropria te theoreti- 
cal description of the electro n scales physics (Schekochi- 
|hin et al. 112009 Bourouaine et al.,_2012) . In the so-called 
TIall-MHD approximation, which is valid in the limit of 
f3i << 1 and Ti < < Te , the Hall term brings a correction 
to the classical MHD theory proportional to di, the ion 
inertial length. This theory can be extended to electron 
scale physi cs within th e inco mpressible Electr on-MHD 
(EMHD) (Biskamp et al. 1999; Meyrand & Gahicr 201(| 
or th e reduced two-tluid-theory (^Sahraoui et al._ 2003 
2012| ). Within the EMHD model, which is appropriate 
to describe whistler mode turbulence, de arises as the 
relevant scale to describe electron scales physics. How- 
ever, in the SW at lAU, where ~ 1 and T./Te > 1 
(see FiglTl), these po pular m odels are not rigorously valid 



( SchekocEihin et al. 2009; Howes 2009), although they 



may reproduce some observed properties of SW turbu 
lence such as the scaling of the magne tic or the electric 
energy spectra in the dispersive range (Matthaeus et al. 
|2008| ). 
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Fig. 10. — Correlation between the frequencies of the observed 
spectral breaks with the Taylor-shifted electron gyroscale pg, 
(right) and inertial length de (left). 

In kinetic theory, 3D recent PIC simulations of electron 
scales whistler turbulence showed that a spectral break 
occurs near dg ( Chang et al.| 2011 1 and not at pe (note 
that in that pa]Der = 0.1, which means that pe = 
del\fW^ ~ 3de)- On the other hand kinetic theory of 
KAW turbulence in the limit fcy << predicts that the 
relevant scales are pi and p^- A simple explanation of 
this can be given by examination of the linear dispersion 
relations of the KAW at different values of /3i (or /?e) 
shown in Fig. 13 (top). One can see that regardless the 
value of Pi the KAW becomes dispersive at kpi ^ 1 and 
no significant change is observed at kdi ~ 1. 

From the point of view of observations answering the 
question as to whether pi or di is the relevant scale is 
more difhcult because, often, SW data at lAU show that 
/3i ~ 1, which means that one cannot distinguish the two 
scales. Nevertheless, so me statistical studies have tried 
to answer t his question ( Leamon et al.|1999 Bourouaine 



et al. 20121. For electron scales only a tew studies based 



on a limited number of events have reported the relevance 
of Pe to charact erize electron scale physics in the SW 
([Sahraoui et al . 2009; Alcxandr ova et al. [2009 Sahraoui 
et al. 201Ua[ ). Alcxandrova ct ain 2012| have conhrmed 
that result over a larger data samples (100 spectra). In 
Fig. [To] we computed the correlation between the fre- 
quencies of the observed spectral breaks with fp^ and 
fd^. The figure shows a relatively high correlation (0.72) 
of the spectral break with the electron gyroscale p^ than 
with the inertial length de - This result contrasts with the 



result in Alcxandrova et al. (2012) which showed no cor- 
relation (C ~ 0.03) betweeii the spectral breaks (given 
by the used double-power-law model) and Pe- A possible 
explanation of the lack of the correlation in that work is 
the downsampling frequency of the STAFF-SA that does 
not allow one to capture properly the spectral breaks. We 
note that the moderate correlation with de (0.58) maybe 
due to all intervals that have /3e ~ 1 (for those events 
Pe ~ c?e)- We note finally that a similar correlation co- 
efficient bet ween /{, and Pe is f ound in magnetosheath 
turbulence (Huang et al. 2013), but a weaker correla- 
tion is found with de ((^' = 0.27). All these observations 
suggest the relevance of pe as a dissipation scale for the 
SW and in the magnetosheath. 

6. DISCUSSION 

The results shown in the previous section underline the 
need to carefully handle the Cluster/STAFF data in the 
SW in order to analyze properly the energy spectra at 
electron scales. We showed in particular that in order to 
evidence the presence of the spectral breaks at electron 



scales time series with high SNR are needed. For the 
selected data a large number of magnetic energy spec- 
tra was found to have a clear break near the electron 
scale, followed by steep power-law-like spectra The 
distribution of their slopes a was shown to range from 
—3.5 to —5.5 with a peak near —4. This distribution 
is wider than that of the slopes in the dispersive range 
which covers the values —3.2 ^ a < —2.3. A compar- 
ison of these distrib utions to those rep orted in magne- 
tosheath turbulence ( Huang et al.[[2013 ) shows some in- 
teresting similarities and ditterences. In the dispersive 
range the distribution of the slopes is very similar: both 
are narrow and show a peak near —2.8, in general agree- 
ment with most of theoretical and numerical predictions 
on SW turbulence. However, in the electron dissipation 
range the magnetosheath spectra were found to be gen- 
erally steeper with slopes a — 7, whose distribution 
peaks near —5.5. This comparison stimulated a further 
analysis of the data to explain the difference between 
the SW and magnetosheath observations based on the 
differences in the SNRs between the two regions. Let us 
examine this effect in detail. 

Although the data used here have the highest found 
SNR in ten years of data survey in the SW, it might be 
possible that even SNR values as high as 10 are not suf- 
ficient enough to determine accurately the scaling above 
fb- The reason is the limited extension of the spec- 
trum above /h (typically less than one decade) due to 
the proximity of the sensitivity floor of the SCM. To 
test this hypothesis we plotted in Fig. [Tl] the measured 
slopes of the spectra above /p^ as a function of the corre- 
sponding SNR. The figure shows a moderate correlation 
(C — 0.53) between the SNR and the slopes. Typically 
most of the slopes \a\ < 3.5, which were not observed 
in magnetosheath turbulence, correspond to SNR< 15 
and the highest values of the slopes \a\ > 4.5 correspond 
to SNR> 15. This correlation suggests that data with 
even higher SNR are needed to fully address the actual 
scaling of the magnetic energy spectra at sub-electron 
scales in the SW. Based on this conclusion the power- 
law fits and the interpretation of the low values of slopes 
(|a| < 3.5) reported in this study should be considered 
with some caution. We emphasize that a weak corre- 
lation (C < 0.1) is found in magne tosheath turbulence 
when SNR> 25 ( [Huang et alT[[2013 |. 

Another interesting point that is worth a discussion is 
the possible dependence of the scaling of the magnetic en- 
ergy spectra in the dispersive and the dissipation ranges 
on the plasma If one assumes that turbulence below 
the ion scale is dominated essentially by high oblique 
KAWs as sugg ested by several SW observations ([Bale et 
'ai 2005; Sahr aoui et al.||2009[ |2010a[ [Salem et al.[|l^01:^[ 
Kiyani et al.|26l3 1 such a dependence could exist. This is 
indeed based on the known result from the linear kinetic 
theory that the damping of high oblique KAWs depends 
strongly on the plasma /3i. This can be seen clearly in 



Fig. 13 (top panel) which shows that the lower is the 



Pi the more damped is the KAW mode (for the same 
angle of propagation). Based on this remark one would 
expect to observe a rapid decay of the energy spectra 
at low Pi, which would translate in to steeper power-law 



or exponentially decaying spectra ( Howes et al . 2011c 
To test this hypothesis we plotted on 1 'ig. |l2| the cor 
relation between the observed slopes below and above 
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Fig. 11. — Correlation between the frequencies of the observed 
spectral breaks with the SNR (see text) at 30Hz of the studied 
magnetic energy spectra. 
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Fig. 12. — Correlation between the frequencies of the observed 
spectral breaks with the plasma fii of all analyzed time intervals. 

fp^ as a function of /3i. The results do not show any 
clear dependence between the values of the slopes and 
f3i. A possible explanation of that is the fact that the 
difference in the kinetic damping of the KAW modes is 
more pronounced between /3i ~ 1 and much lower val- 
ues Pi ~ 0.01. This cannot be seen directly on Fig 
(top). To evidence this fact plotted in Fig. 13 (bottom 
the damping of each KAW mode over one wave period. 
The figure shows that the modes with f3i > 0.5 have 
nearly the same normalized damping rate at all scales, 
which is very different from the dampings of the modes 
at /3i = 0.01 and /3i = 0.1. These latters are shown to 
undergo a strong damping (defined when 2Trj /ujr ~ 1) at 
respectively kpi ~ 3 and kpi ^ 10. This may explain the 
absenc e of correlation between the slopes and Pi reported 
in Fig. [12] since most of the observations reported here 
have l3i>Q.h. Larger samples of SW data would allow 
one to obtain lower values of Pi and to test the scenario 
proposed here. Other explanations are also possible such 
as the presence of non-local effects (due to large scales 
shear flows) that would overcome the strong damping 
at low Pi and sustain the cascade down to the electron 

Another ex- 



(2011c 



scales as suggested in Howes et al. 
planation is that KAW scenario discussed here may be 
irrelevant to interpret the present observations. In any 
case, the origin of the large spread of the slopes ab ove 



reported h ere and in magnetosheath turbulence (Huang 
et al.[2013 ) (as compared to those in the dispersive range 



remains unclear and suggests the lack of universality of 
the turbulence at these small scales. The scaling of the 
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Fig. 13. — Top: KAWs dispersion relations and dampings rates 
normalized to the proton gyrofrequency for = 89.99° and the 
given values of Pi. Bottom: the damping rates of the same modes 
normalized to one wave period. 

spectra at sub-electron scales may depend on local elec- 
tron plasma instabilities. Such dependence was indeed 
found between the slopes of the spectr a at sub-ion scale s 
and the local ion plasma instabilities ( Bale et al.||2009 ). 
Proving the existence of the same dependance at sub- 
electron scale requires having high time resolution mea- 
surements (a few ms) of electron distribution functions. 
Such measurements are not available with the current 
space missions. The NASA/MMS mission is expected to 
meet part of these requirements. 

7. CONCLUSIONS 

We have surveyed more than ten years of the Clus- 
ter burst mode SCM data in the SW and focused on 
the magnetic energy spectra below the ion gyroscale pi. 
We discussed several instrumental caveats that poten- 
tially may influence the estimation of the scaling of the 
spectra at electron scales. We showed in particular that 
high SNRs are needed in order to determine properly the 
scaling of the energy spectra. Once all the caveats were 
considered we found a large number of spectra having 
clear spectral breaks near the Taylor-shifted electron gy- 
roscale fp^ . A double-power-law model has been adopted 
to fit the spectra. The distribution of the slopes below 
the spectral break was found to be narrow and cen- 
tered around ~ —2.8, while that above the /f, was found 
wider and centered near the value ^ —4. The peak of the 
distribution of the slopes below fh is consistent with sim- 
ilar observations of SW and magnetosheath turbulence 



(Alexandrova et al. 2012 Huang et al. 20131 and with 



most of the existing theoretical or numerical predictions 
The spectra above /b are however shallower than those 
reported in magnetosheath turbulence. A possible expla- 
nation of that discrepancy is the difference in the SNR 
between the SW and the magnetosheath. A moderate 
correlation was indeed found between small values of the 
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slopes, typically |a| < 4, with lower values of SNR. The 
reported wide distribution of the slopes combined with 
earlier observations that showed exponential-like scaling 
indicate that the physics of the electron scales is possibly 
not universal and might be controlled by several plasma 
parameters. Based on the assumption that KAW turbu- 
lence dominates at electron scales and on the fact that 
the linear damping of the KAW modes (from the Vlasov- 
Maxwell equations) strongly depends on j3i, an attempt 
to find a correlation between the observed slopes and the 
plasma /3j was carried-out, but the result was not con- 
clusive. 

The Cluster-Search Coil data allow for unprecedented 
in-depth studies of electron scales turbulence in the SW. 
However, as shown in this study, magnetic field data with 
higher SNR is now needed to fully and unambiguously 



characterize the turbulence cascade and dissipation at 
electron scales in the SW. Besides that one needs fur- 
thermore to have accurate measurements of electric field 
fluctuations and high time resolution of ions and elec- 
trons distribution functions, which are not available on 
current space missions. A new mission, TOR, is currently 
being designed to fulfill these requirements. TOR should 
allow in the upcoming years to fully resolve the problem 
of turbulence cascade and dissipation at electron scales 
in the SW. 

This work is part of the project THESOW funded by 
L'Agence Nationale de la Recherche (ANR, France). J. 
De Patoul is funded through the ANR grant. The FGM 

and CIS data come from the CAA (ESA) and AMDA 
(CDPP, IRAP, France). FS thanks O. Gurcan for fruitful 
discussions. 
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